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Section 1. Introduction 



One-loop radiative corrections to the Higgs sector in the Minimal Supersym- 
metric Standard Model (MSSM) have been extensively studied by a number 
of authors [|l], ||, ^, H, |, P] . The motivation for this effort is clear. The Higgs 
sector in the MSSM is strongly constrained (at the tree level, in terms of only 
two free parameters one predicts a large number of masses and couplings of 
the physical scalars) and it may play important role in understanding the 
mechanism of the spontaneous SU{2) x U{1) symmetry breaking as well 
as in discovering (or ruling out) supersymmetry. Indeed, if the soft super- 
symmetry breaking scale is relatively high, say, 0{1 TeV), then the lightest 
supersymmetric Higgs boson may well be the first trace of supersymmetry. It 
is, therefore, important to study in detail the potential experimental signa- 
tures of the MSSM Higgs boson and to understand if one can experimentally 
distinguish the lightest supersymmetric Higgs boson from the Minimal Stan- 
dard Model one, without discovering the other scalar states. 

It is now well established and understood that, in presence of a heavy top 
quark, 1-loop corrections in the MSSM Higgs sector are large. They signifi- 
cantly alter the strategies for the experimental search, even for the top quark 
mass close to the present experimental limit (e.g. the channel A^A'^ 
may be open even for rrit = 90 GeV). In addition, 1-loop corrections intro- 
duce the dependence of the predictions on a number of unknown parameters 
of the MSSM. There have been used three main approaches to calculation 
of 1-loop corrections: 1) the effective potential approach (EPA) |p, 2) the 
renormalization group approach (RGE) |Q (with the heavy sparticles de- 
coupled), 3) diagrammatic calculation of all the relevant Green's functions 

S 1,0- 

The purpose of this paper is to give full details on the introduced earher 
[§, ^] systematic renormalization programme for the Higgs and gauge boson 
sectors of the MSSM. Our approach is a straightforward extension of the on- 
shell renormalization of the Standard Model ||^. The results of the complete 
1-loop corrections to the scalar mass spectrum calculated in this framework 
were given in ref. |^ and those for some of the Higgs boson production cross 
sections at the e~^e~ colliders, with all the 1-loop corrections to the Higgs 
boson couplings included, in ref. 0. Here we also extend and systematize 
those results, always with emphasize on phenomenology of the MSSM Higgs 
production at the e+e~ colliders. Systematic discussion of the Higgs boson 
decays is presented elsewhere 0. 

A systematic renormalization programme for the Higgs sector in the 
MSSM is useful in several ways. Firstly, our calculation provides a "ref- 
erence frame" for the other two methods which use certain approximations 
(e.g. they neglect the contribution of the gauge sector, of the dependent 
effects in the scalar self energies, of 1-loop corrections to the vertices etc.). 



1 



Thus, the accuracy of different methods in various regions of the parame- 
ter space is now well checked. Secondly, our formalism can be consistently 
extended to calculate loop corrections to other sectors of the MSSM. And 
finally, comparison of our renormalization framework with the other two 
methods as well as with the diagrammatic calculation based on different 
renormalization conditions ^ helps to clarify certain points on the choice of 
independent parameters for the Higgs sector in the MSSM and on its gauge 
and renormalization scheme dependence. 

It is well known that the MSSM Higgs sector is conveniently parameter- 
ized in terms of tan/3 = V2/V1 (where Vi are the VEVs of the two Higgs 
doublets present in the model) and the pseudoscalar mass M^o. These two 
parameters are directly related to the free (and unconstrained) parameters 
in the scalar Higgs potential. In addition, tan (3 is an important parameter 
for model building, with theoretical prejudice for 1 < tan/? < mt/mf, [|10| 



However, one should be aware of the fact, that tan/5 is not an observable 
and beyond the tree level it is gauge and renormalization scheme dependent. 
It is, of course, a standard feature of perturbation theory that physical results 
(e.g. predictions for cross sections expressed in terms of measurable quanti- 
ties such as masses and charges) are renormalization scheme dependent up 
to the order higher than considered (i.e. at 1-loop level this dependence is 
a 2-loop effect). However, tan/? (as any other renormalizable parameter in 
the lagrangian), when expressed in terms of the physical masses M^o and 
M/jO (/i° is the lightest scalar), depends on the renormalization scheme at 
the same order at which it is calculated (i.e. at 1-loop level this dependence 
is also 1-loop effect). This implies, in turn, that e.g. 1-loop corrected values 
of M/jO, when expressed in terms of tan/3 and My^o, are (for the same 
values of tan/5 and M^o) scheme dependent to 1-loop effects and the same 
applies to cross sections for the h'^ production and its decay widths. 

In particular, the separation tan (3 < 1 and tan /5 > 1 is renormalization 
scheme dependent already to 1-loop effects and one cannot use the theoretical 
limit tan (3 > 1 unless in the same (or at least effectively similar) scheme as 
the one used to get it. This discussion becomes clearer when we compare the 
calculations in the gauges and in the Landau gauge. As it will be seen 
in detail in the next section, gauges require additional counterterms Svi 
to renormalize tadpole diagrams, whereas in the Landau gauge they can be 
renormalized with the mass and coupling constant counterterms only. The 
additional freedom in renormalization conditions for 6vi can make tan/3 
arbitrarily different (already at 1-loop level) from its definition used in the 
EPA and RGE approaches (which use Landau gauge and are consistently 
employed to get the limit tan/3 > 1). 

Our scheme has been deviced to be very close to the Landau gauge cal- 
culations, although using the 't Hooft-Feynman gauge for easy work with 
Feynman integrals. We can, therefore, use theoretically favoured constraints 
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for tan (3. One should remember, however, that they are based on additional 
assumptions (like demanding radiative gauge symmetry breaking in grand 
unification models) which are much stronger than supersymmetry itself. For 
phenomenological study it is therefore sensible to relax them. On the other 
hand some limits on tan j3 are crucial for a meaningful discussion of loop 
corrections. In particular, for tan/3 0, because of the growing top quark 
Yukawa coupling, the 1-loop corrected M^o oo (eventually perturbation 
theory breaks down). In this paper we study the region 1/2 < tan (3 < mt/mj, 
of the renormalized tan f3 (in our renormalization scheme). The lower limit 
is taken from constraints on the tree level Yukawa couplings [11| - in our 
renormalization scheme the Yukawa vertices receive only generic electroweak 
corrections (large tadpole corrections are absorbed into Higgs boson renor- 
malization constants) and therefore the tree level lower limit for tan (3 can be 
used as a sensible constraint for our renormalized tan (3. The upper limit is 
taken from the results on the mechanism of radiative SU{2) x U{1) symmetry 
breaking studied in the Landau gauge in the effective potential and renor- 
malization group approaches; as mentioned above our renormalized tan (3 is 
effectively very close to the one used in the two other approaches. 

Given the scheme dependence of the {tan (3 , Mj^o) parameterization it 
seems appropriate to complement it with the {M^o, Ma_o) parameterization 
of all the other physical quantities in the Higgs sector. Here there is however 
the well known ambiguity: the tree level mass spectrum of scalars is sym- 
metric under the interchange tan (3 ^ cot (3 whereas their couplings and, in 
consequence, cross sections and decay widths are not. For fixed rrit and m^, 
this symmetry is broken at the 1-loop level by Yukawa couplings which take 
different values in the two regions but the two-fold ambiguity of (M/jO, M40) 
parameterization remains to be present in certain mass ranges. Still, we find 
the {MfiO, M^o) parameterization very convenient for phenomenological pur- 
poses and will use it heavily in this paper. 

The rest of the paper contains the details of our formalism and the last 
chapter is devoted to an extensive phenomenological study for the e~^e~ 
colliders. 



Section 2. Lagrangian and counterterms 

In this Section we define our renormalization scheme for the Higgs sector of 
the MSSM. In general we follow the notation and conventions of ref. 
where the full lagrangian and the complete set of Feynman rules for the 
minimal supersymmetric standard model is given. Let us recall several most 
important for us formulae. We denote the two Higgs doublets as Hi = 
(HlHf), H2 = {HlHl), < Hi >= {vi/V2,0), < H^ >= {0,V2/V2). 
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Scalar kinetic terms read: 

+H2 (d, - '-giB^ - tg2T^W^^ (d, + '-g,B^ + ig2T^W^^ H2 (2.1) 

(T" are the SU(2) generators normalized so that Tt(T'^T^) = ^6"'^). Scalar 
potential has the well known form: 

Vo = m'[H^H^ + ml^2H2-m\^(tabHlHl + c.c) 

+ ^(^7^ + gl)(H,H, - H2H2f + \gl I H1H2 r (2.2) 

where is defined to be negative and ei2 = —621 = — 1. The renormaliza- 
tion constants are defined as follows: 

Vi ^ Z]lf{vi- 5vi), gi ^ ZiZj^^^^gi, g2 ^ Z2Z^^'^ g2, (2.3) 

Hi —>■ Z\lfHi, B^ Zg'^B^, Wfj_ — ^ zly^Wfj, 

Having renormalized all fields and parameters as indicated above, we 
introduce gauge fixing term, expressed already in terms of the renormalized 
quantities (in this respect our approach differs slightly from that employed 
in Hollik's paper |^): 

where 

/ = d^B^ - -^^^g, (hIv2 - Hlv2 + Hlv, - h\v,) 

JT'^ = d^'W'^ - i^W^W^92 {H2T'' < H2> - <H2> T^H2 (2.5) 
+ HiT" < Hi> - <Hi> T'^Hi) 

We choose to work in the gauge " infinitesimaly" different from the 't Hooft 
- Feynman gauge. Therefore we set 

^1,2 (2-6) 
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and next use ' as finite counterterms. As usual, we define physical 
gauge boson fields W^, and by the formulae 

where cw = cos 9w = 92/ \J Qi + 921 sw = sin6'iy. We will also use the 
following combinations of the gauge fixing parameters: 

= cos' 9w + sin^ 9w, a = ^f" sin' 9w + cos' 9w, 



a' - - ) sin 9^ cos 9^, e - (2.8) 

It is convenient to introduce the decomposition of the renormalized Higgs 
fields: 

Hl = -j={v^ + (j)^ + iip,), Hl = ^{v2 + <p2 + ^V2) (2.9) 

and, next, the tree level mass eigenstates for the renormalized fields: scalars 
= (if°, h^), pseudoscalars -f^°+2 = (^°) ^^d charged Higgs bosons 
Hf = {H^, G±), where i=l,2 (G° = and G± = Hf are Goldstone 
bosons). In what follows we will use compact matrix notation introduced in 
[0. The H^'s and -f^f+2's are related to the original fields by the rotations 
|1,[T|]: 

<j^. = ZtHl ^, = Z%Hl, (2.10) 



cos a — sin a \ / sin /3 — cos (3 



where 



V sin a cos a J ' ^ \ cos (3 sin (3 
For the charged Higgs bosons we have 



(2.11) 



For arbitrary values of vi, V2 (not necessarily minimizing the tree level 
potential; they may e.g. minimize the 1-loop potential) the angles a, (3, 9 
are given as follows. Let us first rewrite the scalar potential in terms of the 
fields (pi, (fi and H2, Hf. In this basis the symmetric mass matrices of the 
scalar fields have the form 

L^^^. = 'l <Pi [Ml)^^ 0,- (2.13) 
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and similarly (TW^f)^ for L^^^^ and _ ior Lh2^i. We have: 

(Ml)^^ = ml + \{gl + gl){^vl-vl) 

{Ml)^^ = ml + \{gl + gl){^vl-vl) (2.14) 



{Ml)^^ = + + (2.15) 



The angle a is then given by 

tan 2a - ^ '^'^^^i^ /o 171 

and the tree level scalar masses read: 

;My^^cos^a+ (Al|)^^sin^a+ (A^|)^^sin2a (2.18) 



H 



ml = {Ml)^^ sin^ a + (mI)^^ cos' a - {M^)^^ sin2a (2.19) 

The analogous formulae hold for tan 2/3 and tan 2^. If vi, V2 minimize 
the tree level potential we have — 13, Z± — Zh with 

tan/3 = va/vi (2.20) 

(this is the case with our renormalization conditions). Then m^o = f^%± = 
but one has to add also the gauge fixing terms. Note also that in this case the 
rotations remain unaltered by the inclusion of the gauge fixing contributions 
to the mass matrices. They only shift the Goldstone boson poles to m^o = 
M|, m?Q± — (in the 't Hooft - Feynman gauge). 

The counterterms in the Higgs sector read: 
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a) Counterterms linear in and 



5, Zl E] 



(2.21) 



where 



5i 



v\ dml — ml 6vi + f 2 5m^2 — "^112 

+ 2 Vl V2 Sv2 + (2vf - Vl ^2) SZhi - Vl vl 5Zh2 
and 62 is given by a similar equation with indices 1 and 2 interchanged 
b) Counterterms bilinear in the scalar fields (without derivatives): 



(2.22) 



SLss^-lzi{6Ml) Z^^H^H', 



(2.23) 



with 

{SMl)^^ = 6rni + l{6gl + 6gl){3vf-vl) (2.24) 
+ ^(fl'1+5'2) '2v2Sv2-6viSvi-vI{SZhi+SZh2) + GvISZ^ 



{5M 



+ 



5ml^ - ^ (5gl + 5gfj V1V2 (2.25) 
^ {pi + gt) [v\^V2 + V26V1 - V1V2 {SZhi + SZh^)] 

and (5A^/j)22 given by the equation with indices 1 and 2 interchanged as 
compared to {SMji)^. Similarly: 

SLpp ^-\z'h {^M],)^. Z^^ Hl.Hl, - \ bil Ml ElEl (2.26) 
where 

{^M\)^^ = bm\^\(bg\^bgt)(y\-vt) (2.27) 
+ \ {fx + gt) \^V2?>V2 - 2vi6vi - vl {6Zh, + + '^vldZn, 



{5M'h) 



12 



— (5m?. 
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(2.28) 



and {5Ai']j)22 is given by the equation with indices 1 and 2 interchanged as 
compared to {SAi'^j)]^-^. Finally: 



5L 



± 



Z^^{5Ml) Z^iH-Ht 



(2.29) 
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where 



2 V2 Sv2 - vl {SZhi + SZH2) 



1 2 

1 



(2.30) 



+ o (^2 + 9I) [2v25v2 - 2vi5vi + 2v16Zh, - vj {SZh, + 5Zh,) 
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1 1 
= - Sml2 + ^ Sgl Vl V2- - gl {vi Sv2 + V2 Svi) 

+ V2 {6Zh^ + 6ZH2) 



(2.31) 



and (SAi'^^^ is given by the equation with vi, 5m1, SZh^ and V2, Sml, 
SZh2 interchanged as compared to (6A4'^^^^. In the equations above we 



have used the abbreviations: 



Sgl = gl (2 5Z2 ~ 3 



w 



(2.32) 
(2.33) 



5gl = gl{2SZi-3SZB)^- gfSZB 

(the last equahty holds because 5Zi — 6Zb due to the U{1) Ward identity). 

c) Counterterms to the scalar kinetic terms. Introducing matrix notation 
{5ZH)ij = diag((5Zij^, SZh^) we can write: 



5V = ^Z^^i5ZH),jZi^d,H^d^H', 

+ \ Z'k {bZn\ Z'i d,HUd^Hl^2 
+ Zl{bZn\Zld,E-d>^Et 



(2.34) 



d) Counterterms to the propagators which mix G° with ^ and 

with have the form: 



(Af zl + ^;^) 9^^//^, 



2^gl^gl 



(2.35) 



5Lii. = 7^92^^ Z'^± W^d^Hr + (SCV - S^^) W^d^R- + c.c(2.36) 
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where 



Af = (-l)k^-gl{y^(5Zi-5ZB + 5ZH,)-5Vk) 

- gl{vk{5Z2-5Zw + 5ZHk)-5vk) 



(2.37) 



= (-1)^ 02 {5Z^ - 5Zb - 5Z2 + 5Zw) 



(2.38) 



Af = (-1)^ K (5Z2 - + 5Z^J - 8v^\ 



(2.39) 



Note that when tan/3 = ^2/^1 and 5vi/vi = 5v2/v2 (as will be the case 
in our renormalization scheme) the counterterm for the Aj^ — mixed 
propagator vanishes. 

e) Gauge bosons sector. Counterterms to the gauge boson kinetic terms read 
(^2 = g^'^^'^ - ^^'^'')■. 



SLtn = l{cl.SZn + sl,SZ^)AldlAZ-^6e,Ald^d'^AZ 

+ \ {s'^6Zs + cl.6Z^) ZldlX - \ ^il Zl d'^d^Zl 

+ cw sw {SZw - 5Zb) Ald%Zl - 5ifAld>^d^Zl (2.40) 

+ 5Z^ W^dl^W- - 5CV W+d^d^'W- 

Counterterms to the gauge boson masses read: 
1 



SL±,s = -[{9l + 9l)SX + 2gl{vf + vl)iSZ2-5Z^)]z'^Zl 



+ -gig2(vl + vt){6Z2-5Zw)Z'^Al 



(2.41) 



+ 4^2 



5X + 2 {vl + vfj (5Z2 - 5Zw)] W+^W; 



We have used here the equality 5Zi = 6Zb to simplify the formula and 
introduced the combination 



5X — vf SZffi + vl dZH2 — 2 Vl Svi — 2 V2 Sv2 



(2.42) 



f) Counterterms for the vertices Z Z Hf and A^Z are: 

5Lvvs = \ ^'zzZ'J Z^>^ZlH^ + ^ 9^92 {5Z2 - SZ^) C\ Z'^AIH', (2.43) 



where 



A|z = vk gl {6Z2 - 5Zw) + \ {gl + gt) {vuSZh, - Svu) (2.44) 



and for future use we have introduced the matrix C}^: 

C)j ^ V, Zl + Zl (2.45) 

g) Counterterms for the vertices Z^H^Hjj^2 ^^HfHjj^2 given by: 

= ^^zi'Z^^fAV. Z'.(h'B>^HI 



9i 92 



A% {5Z, - 5Zw) Al (h^.B^H^) (2.46) 



2V^7f + 92 

where = 3^"- 9^, A'lj = Z]^Z^^ - Z^^Z]^ and (a)^^ = diag (Ai, - A; 
with 

A, = gl 5Zh, + 9l (5^2 - + 6Zh,) (2.47) 

Note that ahhough the tree level couplings A^Z^H^ and A^HfHj_^_2 vanish 
in the initial lagrangian, 1-loop corrections to those vertices are divergent 
and their counterterms are given by the above expressions. The apparent 
contradiction with the "theorem" that, if the vertex does not exist at the 
tree level, corrections to it must be finite, can be attributed to the Z° — 7 
mixing. Indeed, as we will show in the next section (eqs. p.4| , |3.6| ), the combi- 
nation SZ2 — 6Zw is entirely determined by Z° — 7 mixed (unrenormalized) 
propagator at zero momentum: 



1 

cwSwM'^ 



5Z2 - SZw = -——^ n^z(O) (2.48) 



which gets contributions only from gauge/Higgs bosons sector and thus has 
the universal character; no new renormalization conditions can be imposed 
on the A^Z^H^ and A^H^H^^2 vertices. 

Section 3. Renormalization conditions and phys- 
ical Higgs boson masses 

Counterterms defined in the preceding section can be used to renormalize 
divergences of the loop corrections in the gauge and Higgs boson sectors of 
MSSM. Their finite parts must be fixed by choosing appropriate renormal- 
ization conditions. The renormalization scheme we are going to define should 
be compatible with supersymmetry. Since we want to calculate measurable 
quantities, we are interested in preserving supersymmetry of the S— matrix 
and not necessarily of the Green functions themselves. With this restriction. 
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our counterterms are sufficient for the renormalization programme provided 
a supersymmetry preserving regularization is used. (We will therefore use 



dimensional reduction which preserves SUSY at least to one loop.) This 
is because the wave function renormalization constants (like SZh^), which 
violate supersymmetry (e.g. the infinite parts of the wave function renormal- 
ization constants of bosons and fermions from the same supermultiplet are 
different when working in the Wess-Zumino gauge with the physical fields 
only), cancel in the 5*— matrix. Below we specify our "on-shell" renormal- 
ization scheme. Quantities with a hat (which are finite) are obtained from 
divergent quantities (without a hat) by adding counterterms defined in sec- 
tion 2. 

1) Gauge bosons sector. 

We write the 1— PI two point functions of the massive gauge boson propaga- 
tors as 

' - ^) (*'^) + ' ^ (P) (3.1) 



and for the photon as 



and require that 



(g^'^e - k^'k") n; [e) (3.2) 



fi;(A;2 = 0) = 

7^e (A;2 = M^) = 

7^e n| (A;^ = M|) = (3.3) 

Ue ti^z {k"^ = O) = 

that is that the tree level and masses are in fact physical ones and 
that the real photon couples without admixture of Z^. Also the photon 
propagator has the residuum equal to unity 0. Explicitly these conditions 
give: 

.z.- ^ ne (n; (o) - 2 co.o„^ + .J"^ '^^i' ) (3.4) 



SZ. ^ T^e I n; (0) . 2 tan...g|p - ^^^l^^ | (3.5) 



iZ, - I n, (0) - — — + -p^^ I (3.6) 



z) 
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6X 



+ 



7^e (M^n; (o) + 2swcw'n.'^z (o) 
K (M^) - l^ni (m; 



^2 _ 2 



(3.7) 



It is also possible to change slightly the finite parts of these renormalization 
constants so that the is renormalized on-shell including the Z° — 7 mixing 



0: 



7^e ni (m|) 



(n^z (M: 



n; (M|) 



(3.8) 



This amounts to small changes in the finite parts of counterterms calculated 
in (0-^. 



The finite counterterms S^i'^ are used to locate poles of the longitudinal 
parts of the gauge boson propagators at M^, M| and 0, respectively: 

1 



w 



-5Z- 



w 



Tie 



(m^) - (m^) + (m, 



(3.9) 



6^1 = — (cos^ OwSZw + sin^ OwSZb 



Tie 



n| (m|) - ni (m|) + ni (m|)] (3.10) 



5^7 = (3.11) 

Note that 6^^'^ should be finite which provides a useful test of the calcu- 
lation. 
2) Tadpoles. 

In order to make contact with the effective potential approach of ref. 
(which uses Landau gauge) we want to work with the fields and which 
have vanishing vacuum expectation values order by order in perturbation 
theory. We therefore require that the sum of the tree and loop tadpoles van- 
ishes. Furthermore, in our on-shell renormalization scheme it is convenient 
to preserve at 1— loop (and higher) the tree level relations between VEVs Vi 
and the renormalized parameters, i.e. we require that the tree level tadpole 
vanishes by itself. Equivalently, tadpole counterterms cancel the 1-loop tad- 
poles. In Lorentz gauges, where ^2'^ = (of which Landau gauge with, in 
addition, = is only one example), one usually takes 6vi = 0. It would 
be therefore convenient to use Svi to cancel only this contribution to tadpoles 
which appear when ^^'^'s are switched on. However such separation turns 
out to be impossible due to the spurious infrared divergences which such a 
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procedure introduces. Nevertheless, it is possible to do that for infinite parts 
only. We fix therefore 5vi (in the dimensional reduction we employ) by the 
formula: 



5vi 



1 



4 {AtiY 



d-A 



7 — log 477 



(3.12) 



and cancel the remaining parts of the tadpoles with the counterterms 5m\ 
and 5m\. This means that we eliminate these counterterms everywhere in 
favour of tadpoles of the original scalar fields (pi. 



(3.13) 



Notice, that we have bvxjvx = 6V2/V2 = 6v/v, which is essential for the 
definition of the electric charge by the Thomson limit. Without this equality 
the on-shell photon would mix with and the Ward identity ^ leading 
to the correct definition of e would be spoiled. 

After some rearrangements and substitutions, Smf and 6ml can be writ- 
ten down as: 



5m\ 



tan /3 bvr?Y2 + (^1 + 



m 



6v 

V 
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6v 

V 



9i 



cot p om^2 + ("^2 + cot p m^2) — 



(3.14) 



g2 {26 Z2 - 36 Zw) - gl6ZB 
\ [9I + 9t) (vl + ^vl - 2 tan^ {iv\ 



^2-^1 



6v 

V 



+ vl (cot^ /5 - 3) 6Zh, + (2 - cot^ p) 6X 

3) Mixing of pseudoscalars and G° with gauge bosons. 

Because the mass of A^ is usually chosen to parameterize the Higgs sector 
of the model, we require that the on-shell A^ does not mix with Z^. Writing 
the corresponding mixed propagator as [| 



— T)^ V 

P ^ZP 



1,2 



(3.15) 



and demanding TZe S^p (M^ 
6Zhi 



we get: 
,2 



1 , 2 n Jvi 2 Jv2 gi 6X 
1 + cos /3) — + sm 13 — + -TTrT 
Vl V2 4 



tan /? 



■Re Y} 



ZP 



(Ml) (3.16) 



^This is as a compactified version of the standard notation: S^p ^ '^z"A"i ^ip 



^12 
-'SS 



S^ojyo etc. 
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6Zh2 = COS p— + (1 + sm (3)— + -jj^ - J^^^ ^zp{Ma) (3.17) 

The finite " counterterms" can be used to cancel the j — mixing at 
p2 = 0: 

Tie S^p(O) - swCwMz {6Z2 - 6Zw) - (d^f - d^f) = (3.18) 

4) Pseudoscalar propagators. 
To determine 5m^2 '^^ require: 

7^e t]^p (m1) = (3.19) 

which gives: 

(5mi9 = sin 13 cos (3 [Ue S^V {m%>) - Ml {sm^ I35Zhi + cos^ (35 ZH2j (3.20) 



+ sin (3 y- [g^ + g2) {vidvi - V26V2) + {m^6vi + m^2^V2 - T^^ j /vi 

+ cos^ (3 Q [pi + g"^ {V2SV2 - vi6vi) + [ml6v2 + ml^Svi - T^^ /v2 

In addition we can arrange for 

Tie tfp (m|) = 0, Tie E^^_ (m^) = (3.21) 



This, together with ( |3.18| ), determines the three independent constants 6C,2, 
Sir and 5^2^. 

5) Renormalization of the standard model is usually performed in such a way 



[§, that e = (71(72/ y gl + gi is equal to the electric charge as measured in 
the Thompson limit, that is e^/47r = 1/137.036. This is achieved as follows. 
Denote by Zieft and Zright the wave function renormalization constants of the 
left and right handed electrons, respectively. Zi^ft and Zright can be used to 
give the residues of the poles of the left and right handed electron propagators 
equal to unity. Then, writing the corrected electron-electron-photon vertex 
as 

ie-f'' + ieA^ = ie7^ + ieA'^y + ie-f^A'^A (3-22) 

one has 

A'^y = A^^y + ^{6Zi,ft + 6Z„gM) + \{6Z2-6Zw) 



1 1 n'^ Co) 

Af, + -(.Z,„ + ^Z„,„)-— ^ (3.23) 
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where we have used ( p. 4] , ^.61 ). This is the equation which, in the Standard 
Model, ensures that g. 111: 



A':v[q=0,pi=pi = mi)=0 (3.24) 

and the same condition can be used in MSSM, as a result of the residual 
f/(l) gauge symmetry. This completes our renormalization scheme. 

Basic formulae used to calculate the unrenormalized two point Green 
functions, which are necessary to determine the counterterms are collected 
in the Appendix A. Using expressions displayed there and following the pre- 
scriptions given above one can compute all the renormalization constants in 
the gauge and Higgs sectors. They can subsequently be used to renormalize 
all relevant self-energies for arbitrary momenta p^. The cancellation of diver- 
gences between counterterms and bare two point Green functions provides a 
very strict test of the calculations as well as of the computer programmes. 

Having renormalized self energies one usually calculates one loop cor- 
rected masses from the "perturbative" formula: 

Mlo = mlo + Tie t% (p^ = mlo) (3.25) 

i i \ i / 

However, since in MSSM the corrections to Higgs boson masses are expected 
to be large, we do not use this formula and instead determine the masses as 
the exact poles of the scalar Higgs boson 2x2 matrix propagator [Q. We 
thus have M^o = where the two solutions of the equation: 



Ue 



(3.26) 



In fact, we have checked numerically that the 1-loop masses computed from 
the "perturbative" formula differ significantly from the exact poles of the 
propagator. The difference between the two is particularly big for small 
tan/3 as well as for large tan/5 and ~ M^. Similar difference has been 
also noticed recently in ref. 0. The results for scalar Higgs boson masses 
obtained within the approach just described have been extensively discussed 
in ref. 



Equation ( 3.26 ) can be solved only numerically. However it is possible to 
give approximate explicit expressions for the physical Higgs boson masses. 
Their accuracy, as compared to the exact results, is typically 1-2 GeV. They 
read: 

Mlo ^ — r 3.27 

^1.2 2(l-2£) ^ ^ 

where 

a = V(l - 2^) ('^i - '^2)' + {e {61 + ^2) + 263 (1 - e)) 



2 

b = {1 - e) {61 + 62) + 2e63 (3.28) 



15 



and 



Si = Tie T}^s (0) + m]jo 



62 = ne ±fs{0)+mlo 
(0) 



Tie 



(3.29) 



e = 



w 




(3.30) 



e determines the dominant correction to the p"^ dependence of the scalar 
self-energy. 

The following important point is worth emphasizing. Defining the masses 
as the exact poles of the full 1-loop 2x2 matrix propagator means that a 
summation has been performed over all 1-P reducible self-energy diagrams 
with 1-loop segments. This accounts for the leading higher order contri- 
butions to the Higgs boson mass corrections (e.g. in the second order this 
contribution is 0{Y^^) whereas the irreducible 2-loop corrections are at most 
0{Y^) where Yt is the top quark Yukawa coupling). Another remark is that, 
in our renormalization scheme, the on-shell residues of the Higgs scalar prop- 
agators are not equal to unity. Thus, finite renormalization constants for the 
external particles will appear in the S'-matrix (see next Section). Again, they 
are calculated from the 2x2 matrix propagator, i.e. with the 1-P reducible 
diagrams included. 

Section 4. Cross sections for the Higgs par- 
ticles production in the e+e~ colli- 
sions 

The main sources of the supersymmetric Higgs particles in the e'^e^ collid- 
ers are the Bj0rken process e~^e~ Z^h^{H^) (e+e^ Z^Hf in the ma- 
trix notation) and the associated production of the scalar and pseudoscalar 
e+e~ A°h°{H^) (i.e. e+e" 0- The set of diagrams we take 

into account for the Bj0rken process e^e~ Z^H^ is shown in Figure 1 
(boxes represent loops of all possible virtual particles). Exactly analogous 
set has been taken for the second process e+e~ — ^ H^Hj_^2- the present 
level of accuracy we are interested in including radiative corrections only in 
the Higgs sector and we neglect genuinely electroweak corrections such as 
corrections to the electron self-energy, to the Z^e^e~ vertex, and the box 

^ We do not discuss here the possibiUty of production of the Higgs particles by 
bremsstrahlung off heavy quarks (e.g. e~^e~ — > hbhP{H^)), which can be significant for 
large tan (3 . 
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diagrams. Further simplification can be achieved when calculating matrix 
elements: contribution from the diagrams le, f, i and from the longitudi- 
nal parts of the and 7 — Z° mixing self energies in diagrams Ic, d, g 
are proportional to the electron mass (or vanish). Therefore, in the (very 
good) approximation mg = 0, the numerically important corrections can be 
summarized as follows. 

1) . Corrections to the Z^, 7 and 7 — Z° self-energies on internal and external 
lines. They are usually small, of the order of few percent. The only exception 
is the vicinity of the Z° resonance pole. 

2) . Corrections to the scalar and pseudoscalar propagators. They are the 
main source of the differences between the tree and 1-loop results. As dis- 
cussed in Section 2, they are responsible for large changes in the physical 
masses M/jO and M^o which are predicted for given values of tan/? and Ma, 
and also for finite renormalization constants Z^^^ for external particles in the 
S'-matrix elements. The latter are very important as they constitute the link 
between the angle a which enters the tree level Higgs boson vertices in the 
EPA method and our tree level angle a. 

The Z'^^^ factors are given as follows. Let i' denotes the index of the 
"supplementary" Higgs particle: i' = 2(1) when i = 1(2) (formally i' = 3 — i). 
Then the finite external line factor reads: 



Tie 



1 - 



_d_ 



\ 



m 



HO. 



^SS 



(g2 



(4.1) 



and in the case of no mixing on the outgoing Higgs boson line the S'-matrix el- 
ement is given by the truncated Green's function multiplied by the (^|f*)^''^. 
In presence of mixing i' ^ i on the external line (with the scalar i being the 
final state particle) the respective factor on the external line reads: 



A/2, 



(4.2) 



Therefore, the truncated Green's function is in this case effectively multiplied 
by 



1/2 



y Si J Smixi ^ ' ^ 



In these expressions masses denoted by small letters are the tree level mass 
parameters and capital letters stand for the physical masses. Quantities with 
hat usual, renormalized quantities. 
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To compactify the notation it will be useful to treat formally the mixing 
constant as a two-element array: 

= ) (4.4) 

iJrrnx'' \ ' 'J mix'' J ^ ' 

Exactly analogous equations hold for the pseudoscalar constants Z|>^* and 

^, but in that case their effects are numerically much less important. 
3). Vertex corrections. For some choices of parameters they can be numeri- 
cally important, although only O ((^^m^/M^). At the tree level the relevant 
Higgs boson vertices read (assignment of momenta is shown in Fig.l): 

= ^^^'L (P - ^ ip - qY (4.5) 



After including 1-loop corrections and adding the counterterms defined in 
Section 3 the renormalized vertices have the form: 



^ZZS 



'g'" (y^fs + ^i) + P'P'F^ + q^q'F^ + p^q^F^ + q^p''F[\ (4.8) 



where Fa are vertex formfactors. 

The vertices V^'^^ are obtained from the tree level vertices V^^^ by inclusion 
of the 1-loop corrections on the external lines and read: 

T>(0)ii _ / ■p'ext p-exA 
^ZPS — y^^Si ^Pj ) ^ 

Y ZPS Smix^ ZPS Pmixj ZPS Smix'^ ^mixj ZPS J 



(zifZ^fY^' V Zf* ,Jp^\ ^^'^^^ (4.9) 

a,/3=l 



where 



i(a) = |' . "J (4.10) 



and similarly for j (/3), 

1/2 
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with 



Zf^* = 7^e 1 



_d_ 



^1, 



(4.12) 



To include the photon exchange in the s-channel also the knowledge of the 
vertices V^ps and V^zs (vanishing at the tree level) is necessary: 



(4.13) 



g^^Gl + p^p^Gl + q^q^Gl + p^q^Gl + q'^p'^Gl 



(4.14) 



The 1-loop contributions to the vertices Vzps, Vzzs, V-yPS and Vyzs can be 
found in Appendix B. 

The cross sections (in the CMS system) for both processes have the form: 



da 



[Ai + Ai cos^ QcMs) 



A {s^ ml, — ^Js'^ + mf + 1712 — Ssmf — 2sm2 — 2mfm2 (4-15) 



where 



Ai + A2 cos^ 0CMS = ^ E (MM^) 



(4.16) 



M% = evihh, [V,% (cv - CAl') + + V§fZ,z] u{h) (4.17) 



M^s = ev{h)Y [Vz'zs {cv - ca^') + Kzs + V^Ts^'Z^z] u{h)e,{p) (4.18) 



is) 



(i-fi; 



(4.19) 



In the above expressions u{ki) and v{k2) are spinors of the incoming electron- 
positron pair, e^(p) is the polarization vector of the outgoing Z^, Cv = 
(— 1 + 4s^) /AswCw and ca — —1/AswCw- ^'ps and V^zs related to 
V^pg and as follows: 



^iPS 



^iPS 



s-Ml- ni(.) 



(1 



n;(^ 



= pi'G'l - q^G'i 



(4.20) 
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s-Ml-Iil{s) 



s ( 1 - n; (s; 



g^'^'Gl + ^^^^2 + qy^'Gl + p^g'^Q + q^p'^Gi (4.21) 



Performing the sum in eq. ( f4.16| ) one gets for the associated scalar and 
pseudoscalar production: 



1 



''PS 



+ 



4s 



where 



ttpg 

4s 



(o)^i + F^j + FiJ) + s^i"?^^' Z^z + G% + G% 



OA {2%fi +Fii + F^ 



Analogous expressions for the Bj0rken process are more complicated: 



A' 



2ZS 



32M| 



X 



+ 



2a^zs 



s + Ml- Mlo) h- 

i / 

2dzs -(s + Ml- Ml.) 4 



-4sM| 



h'zs 



+ 



-4sM| 



"-zs 



where 



•^izs 




^^zs = 


cv 


'^zs — 


cv 


'^zs = 


ca 


r/* — 


ca 







lO 




lO 


( 




+ 


^zs 





A 



2ZS 



\^^s+Fi)+vm 



ZZS ^iZ + t^i 
4 



(4.22) 



(4.23) 



(4.24) 
(4.25) 



(4.26) 



As it has already been discussed several times, the expressions for the 
cross sections displayed above are not exactly 1-loop formulae: they have 
been designed to include the leading terms in the higher order corrections. 

Section 5. Results 



In this section we recall and extend the results obtained in ref. |^. We 
focus on the predictions for the mass spectrum for the two neutral scalars 
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and and on their productions rates in the e^e~ colhders. Two related 
questions are discussed: the detectabihty of the scalars and as a func- 
tion of the CMS energy and the possibihty of distinguishing between the 
supersymmetric and the minimal standard model (SM) Higgs bosons on the 
basis of the production rates. 

In Fig. 2a the allowed mass regions in the (M^, M^) and (M4, M^) planes 
are shown for rrit = 120 and 180 GeV, and for several sets of values for the 
soft supersymmetry breaking parameters. For the chosen set of the central 
values of the parameters we also show the contours of fixed tan/5 (Fig. 2b). 
The soft supersymmetry breaking parameters are: universal soft masses for 
the three generations of squarks and sleptons, M^g and M^i, respectively, 
gaugino mass Mgau (with the SU{2) and U{1) gaugino masses related by 
the grand unification formulae Mg^^^ = ^ tan^ OwM^^^'^^), universal dimen- 
sionful parameter A in the trilinear scalar couplings AYj (where Yj is the 
Yukawa coupling for the / generation of fermions; f = u,d, I) and the Higgs 
mixing parameter /i in the superpotential. The general features of the results 
shown in Fig. 2 are already well known. Here we illustrate their dependence 
on the values of the above set of parameters and on the top quark mass. An 
interesting observation is that changing fi and A affects the mass spectrum 
differently for different values of the top quark mass and tan/3. Larger values 
of fi and A can increase or decrease the corrections to Mh and Mh, depend- 
ing on the values of rrit and tan/?, what can be understood as the effect of 
squark mixing. The difference between our results and the results obtained 
by the EPA or RGE is typically (9(5 GeV). A more detailed comparison has 
been presented elsewhere 0. 

Two more remarks on Fig. 2 are important. One can see the well known 
2-fold ambiguity: in part of the parameter space the same pair of values 
{Ma, Mh) corresponds to two different values of tan/?. Thus, in two different 
versions of the model with different values of tan/5, scalars with different 
couplings (i.e. different cross sections, decay widths) happen to have the 
same mass. This has the following implication: as seen in Fig. 3, plotting the 
production cross section for the as a function of Mh we get two branches 
of the predictions, depending on the value of tan/?. 

Secondly, there is also another approximate degeneracy: for large Ma 
(> 100 GeV) and large tan/5 (> 5) Mh is almost constant. As seen in Fig. 3, 
in this region the production cross section in the process e^e~ — >■ Z^h^ also 
becomes independent of Ma and tan /?, and to a very good accuracy is equal 
to the cross section for the production of the SM Higgs boson with the same 
mass Mh- Similar remarks apply to the region tan/? < 1: for fixed value of 
tan/?, the mass Mh and the cross sections for the e~^e~ Z^h^ almost do 
not depend on the value of Ma- Also, for large Ma, the dependence on Ma 
of the Mh and of the cross sections is very weak for any fixed value of tan (3 
(and in addition, as said above, the dependence on tan j3 disappears for large 
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tan/3). Finally, for large tan/? (> 5) and small values of Ma (< 100 GeV) 
the mass of the heavier scalar is almost constant. All those features are 
important for experimental analysis: on the one hand the variables Ma and 
Mh (or Mh) are more direct and physical to use (in particular, one should 
remember about renormalization scheme dependence of tan j3). On the other 
hand they may be inconvenient in certain regions of the parameter space and 
the use of the unconstrained parameters (M4,tan/5) provides a useful "blow 
up" of those regions. The complementarity of both parameterizations will 
be further illustrated below. 

In Fig. 3 we show typical behaviour of the cross sections for the processes 
e^e~ Z^h^ and e^e~ — >■ A^h^ as a function of the M/i, for fixed values of 
Ma- The values of tan/3 change along the curves, as marked on the Figures. 
We see the well known complementarity of the two cross sections, which 
follows from the tree level couplings and is preserved by radiative corrections 
(as is clear e.g. in the EPA): the pair production is dominant for small values 
of Ma and large values of tan/3 , i.e. in the region where the cross section 
for the Bj0rken process is very small. 

In Fig. 4 we show the regions in the (M^, Mh) and (M4, Mh) variables 
in which at least one of the cross sections for the two production processes 
considered here {e^e~ Z^h^ and e^e~ A^h^) is larger than do, for 
different values of ctq, different values of the CMS energy and the top quark 
masses = 140 and 180 GeV. In Fig. 5 we show analogous plots in the 
(M^,tan/3) planes. It should be pointed out that for -^i = 500 GeV there 
is one additional very important production process, namely the WW fusion 



1^, which we do not consider here. The experimental signatures for the 
Higgs boson production via this mechanism are, however, different from the 
signatures in the two processes discussed here, so one can consider them 
independently. 

The effective regions in Fig. 5 are superpositions of the regions where 
(^Ah{H) or azh{H) are dominant. For instance, in the (M4,tan/5) plane and at 
high enough energies, the observed structure with the dip for Ma = 100 GeV 
reflects the complementarity of the two cross sections: the first one dominates 
to the left of the dip (apart from the very small values of tan/?) and the 
second one to the right of the dip. Notice that for realistic values of ctq the 
region of unobservability of the h*^ does not disappear with increasing energy. 
Fortunately, in this region of the parameter space the heavier neutral MSSM 
scalar becomes observable at high enough energy. For rrit = 140 GeV, 
max{azh, o'Ah) or max{azH, cfah) are larger then 0.05 pb in the whole range 
of (M^,tan/5) for y/s = 200 GeV and for rrit = 180 GeV the same occurs 
for y/s = 240 GeV (for a broad range of values of the soft supersymmetry 
breaking parameters). It is also worth observing that a small increase in 
energy can radically change the range of {Ma, tan jS) which can be explored. 
For instance, for rrit = 140 GeV, the increase from ^/s = 190 GeV to ^/s = 
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200 GeV gives the result shown in Fig. 5a. This is easily understandable as 
due to the discussed earlier approximate independence of the Mh on the Ma 
and tan (3 for large values of those variables: changing tan /3 from, say, 4 to 
50 results in only a tiny increase in the value of and a small change in 
energy is sufficient to reach the discovery limit. 

Finally, in Figs. 4 and 5 we show also the regions where the cross sections 
(^zh{H) are very close to the cross section for the SM Higgs boson production. 
This actually occurs for most of the yet unexplored parameter space. Similar 
conclusions hold for the decay branching ratios [|l8l , |^. This is the general 
property of the MSSM: for a sufficiently heavy pseudoscalar A, the effective 
low energy theory (after decoupling of the heavier Higgs scalars) contains 
only one SM-like Higgs boson. Detailed inspection of Fig. 5 shows that 
for mass ranges not excluded yet by the LEP measurements and for the 
values of (Tq considered here we have two possibilities in general. First, only 
one MSSM neutral scalar is observed via in the Bj0rken process e~^e~ — >■ 
Z^h^{H^) and is almost always indistinguishable from the SM Higgs boson 
(at least on the basis of the cross section measurements). Second, both the 
lighter and the heavier scalars are observed, giving us strong direct evidence 
for supersymmetry. 

The authors wish to thank A. Dabelstein for careful checking of the large 
part of the formulae and W. Hollik for helpful discussions. 

The formulae for self-energies and vertex formfactors presented in the 
paper have been prepared as the library of FORTRAN procedures. They are 
available on requests at the e-mail addresses: 
rosiek@fuw.edu.pl (or rosiek@plearn.bitnet) 
chank@padova.infn.it (or chank@ipdinfn.bitnet) 

Appendix A. Two point Green's functions 

In this Appendix we display the formulae for the unrenormalized two point 
Green functions. In order to preserve supersymmetry of the S*— matrix ele- 
ments all Green functions have been regularized using dimensional reduction 

[Hi- 
First, we define our conventions for vertices which appear in the formulae 

below. We write vector boson(/i,a)-incoming fermion(l)-outgoing fermion(m) 

vertex as 

- i -f" (4"^ - 7^ cf^) (A.l) 
For scalar(i) replacing vector boson we write 
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where Pl,r are left /right projectors 0. Interactions of three and four scalar 
fields in our convention take the form (eg. scalar two pseudoscalars vertex 
and two scalars two down-type squarks vertex): 



• -trijk 
^ ^ SPP^ 



^ ^SSDD 



(A.3) 



and couplings of vector boson to sfermions, e.g.: 



• -trmn 

* ^VDD 



(k + pY 



+ i vr"" 



V1V2UU 



(A.4) 



(where k^^ is the incoming momentum of squark). The specific form of 
all vertices can be found in ref. [|12[ Q 

Our notation for (s)particle masses is as follows. We denote lepton, d- 
and u- type quark masses as mi^, m^^ and respectively where / is the 
generation index (/ = 1, 2, 3). Similar index I is used to count 3 generations 
of sneutrinos whose masses are denoted as rricj. Slepton and squark masses 
are written as m^^, mD„ and mu„ and indices n vary from 1 to 6. 
The two chargino masses are denoted as mc\ {i = 1,2) and the neutralino 
masses as rriN^ (j = 1,2,3,4). Let us stress that all masses we use are 
the true eigenvalues of the corresponding mass matrices. Rotations from the 
initial fields to the mass eigenstates are described in |]12| . 

We define also some auxiliary functions which help us to make the for- 
mulae for self energies more concise: 



B{p^, ml, ml 



ao{ml) 



1 2\ ^1 ~ ^2 

- an(ma) — 

3 °^ 3p2 



ao{ml) 



ao{ml] 



+ ^o(p^ ml, ml) + ^(^ml + ml-^ 

(p^ — ml + ml)^ 



3 p^ 



&o(p^ ml, ml) 



(A.5) 



For scalar contribution to Tl^ (see |A.14| ) we define S^{p 



ml,m2 



S 



B 



4 B{p^ , ml, ml) — 2 ao(m^) + 2 aQ{ml) 



+ [p^ — 2 ml + 2 m^ &o(p^; ml, ml) 



■^Notice the difference with the ref. 



where this vertex is denoted as 



(A.6) 



4 



Notice that expressions for neutrahno couphngs to scalars and pseudoscalars given 
in ref. ||l2| should be symmetrized (without dividing by 2) in fermionic indices. Also the 
factor 2 multiplying 6^^ in the coupling to charginos should be removed and in the 
neutralinos vertex factors multiplying left and right projectors in the bracket should 
read: Zfj*Z^J 



and Z-'^Z 



3i 



Z^Z"^* , respectively. In expressions for the 



sfermion mass matrices one should change the sign of the terms proportional to 
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and for fermionic contribution to 11-^ and 11'^ we use F'^^{p'^,m'f,m'^): 



pab 



-2 {cfc 



v 



') [4622 (^»^ , m'^j) + aoim'^i) + ao("^j) 



(A.7) 



B 



A ( aij hji I aij bji 



2B{p\mlm])-a,{ml) + a,{m]) 



+ [p^ - + mfj bo{p'^,m'^,'m^) 



(A.8) 



Chargino and neutralino contributions to the Z°-pseudoscalar mixing can be 
expressed in terms of 



7^/ 
Ci 



(\( alk I alk\ ikl ( alk alk\\ ikl 

— \[[cv +Ca )Cl -[cv -Ca ))Cr 



p^ 



X 

+ 

X 



ao(mf ) - ao{ml) - (p^ - mf + ml) boip"^, ml, mf) 



fllk „a,lk\ ikl 



V 



(alk 1 alk\ ikl 
\Cy + Ca j Cp, 



rrik 



ao{mf)-ao{ml) + {p'^-ml + ml) 6o(p^ "^L "^f)]) (A.9) 



Finally, for chargino and neutralino contributions to the scalar and pseu- 
doscalar self energies the auxiliary function is useful: 



pff 



^ ( ikljlk I ikljlk\ \( 2 2 2\ j / 2 2 2\ , I 2\ 

^ \ Cl c'r +Crc'l ) [[V -m^jboip ,mi, ^k) + aoi^k) 



+ ao(mf )] - (cf cL'' + 4'4') m b,{p\ , ml 



(A.IO) 



In all expressions here and in the other appendices we use the standard 
Veltman functions oq, bij and Cij defined as follows 0: 



d<^k 



{2tiY - m2 (47r)^ 



aQ[m 



(A.ll) 



i-d 



1, fc^, kf'k^ 



{2ttY [k^-ml] [{k + pf-m 

bo{p^,ml,ml) 
piJ- hi{p'^, mf, m\ 



I ) I pi^p" h2i{p'^-,m\,m\) + gf^" b22{p^-,m\,m\) 



(A.12) 



^We use the standard metric (H ). 
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{k+pY 



rrir 



{k +P + qY 



(A.13) 



—I 
(4^ 



' Co{p,q,ml,ml,ml) 
P^ cii(p, g, mf , m|, m|) + Ci2(p, Q', mf , m|, m|) 
P^p" C2i(p, g, mf, mi, m|) + g^g^ C22(p, g, mf, mi, m^) 

+ {pf^q" +p''q'') C23{p,q,ml,ml,ml) 

+ g'"' C24(p,g,mf,mi,m|) 



1) Vector bosons self-energies (and propagators which mix gauge bosons). 
Using the decomposition 



(A. 14) 



(where U^{p'^) = [n^(p2) - n^(p2)] /p2 ^j^g notation of Section 3) we 
have: 

a) Z° boson self energy (we suppress the momentum argument in all func- 
tions) : 



(47r)^ni = - 14622(0,0) +p' 6o(0,0) 
+ E {ml , m^ + 3 Fr {ml , m^ J + 3 Ff^ {ml, ml 



Ul 1 Ul 



1=1 



m,n=l 



1=1 



+ 3 iV^DDr &22(m^^, m^J + 3 iV^JVr b22{ml^,mlj^ 
+ E ^iz.. Mml) + E [^zzLL «o(miJ 



7=1 



n=l 



+ ^VS^^j,ao{mlJ + 3V^^uu(^o{ml^ 



(A.15) 



+ E O + ^ E Fr{m%,m%; 

2 

+ 2e^ cot^ 26^14/ E 1 2 ^22 ("^i+ , m't+ ) oq {m't-^ 



+ 



2 2 



+ ;;f E (^m) b22{mlo,mlo ] 



2 4 



ij=l 

.2„2 nj2 



_92 

4c^ ,=1 



Eao(m^o) 



.4 2 



-W i=l 



2glcl, [2al{Ml) + (2^^ + M^) 6o(M^, M^) + 4622(M^, M^) 
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(47r) = 

+ E [Fi' (mlml) + 3 (m^^, m^J + 3 J 



45(0,0) + 6o(0,0) 



7=1 

3 



+ (ml, ml) + E 

' ^ m,n=l 



+ 



+ 



a,n=l 

' ^ ' 5\VzTu\'S^{ml^,mlj] (A.16) 



+ e 



2 2 ^ 



b) Photon self energy. In this case we must have p^II:^ = — so we write 
only n^: 



(4^)' = E [ft [K^ r^A) + 3FX' (m^^, m^J + WT ( 



7=1 



6 r 2 4 

+ 46^ E h2i{m\^,m\^ + - 622(m|)„, m^„) + - 622 m^J 

n=l L 
6 



n=l 



(A.17) 



+ E (^c. , ^cj + 2e' E [2 ?>22 , ) + oq (m^+ ) 

i=\ 1=1 

+ 4e2 [ao(M^) + (p' + 2M^) 6o(M^, M^) + 2 622(M^, M^) 
c) mixing with the photon: 

(4vr)^ n^^z = E O + 3i^/ K, m^J + (m^^, m^J 

7=1 

- 4e E [v-i^r^ 622(mi„,miJ + V^J^^d ^'22(m^„,m^J 



n=l 



- 2 V1'«^622«„,m^^J 



(A.18) 



+ E ^^ZLL ao(miJ + 3 ao(m^J + 3 V:;^^^ ao{mlJ 



n=l 
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i=l 

+ 2eg2Cw [2ao(M^) + {2p^ + M^) bo{M^, M^) + 4622(M^, M, 

2 

+ 2e^ cot 26l|^ V [2 622 , ) + «o ( 



■2 ^ 



[m 



{^^f n,^z = E O + 3i^/ K, + 3F/ {ml.ml) 

1=1 

n=l 

2 yj^^^ S^iml^,mlj]+j:Ff{ml,ml) 



(A.19) 



+ cot 2ew E 'S'^("^^+, "^^+) - e 5^2 cwiri \^ «o(M^) 



i=l 



+ (l0p2 + 8M^)6o(M^,M^) + 8M^--p^ 
d) boson self energy : 

(4^fn^ = E[Fr^K,0) + 3F7^(ml^,m: 



3 6 



+ ^j:j:\VwoL\'b22iml,mlJ + i'^ E I^t^wI' &22(m^„, m^^J 

7=1 n=l n,m=l 
3 6 

+ E ao(m?J + [VwwLL ao(miJ + 3 V^V^^j oo(m^J 

7=1 n=l 

+ 3 y^Vf/c/ «o(m^j] + E E W,, 

i=i j=i 

2 .. 2 ^ 
+ ^2 E (^m) ^22 (m^o , ) + ^2 E ^22 (m^cY^ , ) 



(A.20) 



i=l 



1 ^ 1 2 

+ 7 ^2 E «o(mHo) + - ^2' E «o(m^+) - M^. 60 (M^,, 0) 
4 i=i '2 



^ 1=1 

+ ^2 [2ao(M|) - ao(M^) + (4 + m| + M^) bo{M^, M|) 
+ 8 622 M^)] + 3 ^2' ao(M^) 



(4 p2 + M^) bo{M^, 0) - ao(M^) + 8 622 (M^, 0) 
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K = i: [K"^ H, 0) + 3 Fr «, <) 

1=1 



3 6 



1=1 n=l n,m=l 



=1 i=i 



1 1 r 

+ E ^^+) - ^2 ^ [4 ao(M^) + 4 ao(M|) 



- Ml 



a^{Ml) - ao{Ml)) - 8 



(M, 



w 



m: 



2^2 



+ (l0p2 + + 2M|) 6o(M2„ M|) + + 4M| - ^ p- 

o 1 



6o(M^,M|) 
4 



— e 



3p2 



/l/f2 /M-4 

4 ao(M^) - 8 ao(M^) - 8 bo{M^, 0) 



+ (lO/ + 2M^) 6o(M^,0)+4M^--p2 



(A.21) 



2) Mixing of vector bosons with pseudoscalars. We denote the Green's 
function of the incoming Z° with four-momentum p'^ and the outgoing 
pseudoscalar as: 



Using this convention we have: 



(A.22) 



3 ^ bo{ml ml ) 



- 3 Z^jj ^ bo{ml^,ml^) 

V2 



+ i E K-L^/rr h(mij - ao(miJ + K - mi J &o(mL, mi J 

" ri.m=l 



I Q T/nm Trjmn 



MmlJ - aoimlj + (ml^ - ml J bo(ml^,mlj]) 



8 E C}, At 1 [(3 - M| + m^o) 6o(M|, m^o) 
- ao(M|) + ao(m^o)] + E Vf^^ 1 k(m^o) 
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i,fe=i 



l,k=l 



(A.23) 



The matrix Z]j is defined in Section 2. 

In the same convention, for the 7 — G° mixing we get: 



(A.24) 



At 1-loop there are no graphs which mix 7 with A^. Thus, if the 
finite part of the counterterm to this mixing vanishes (as is the case in our 
renormahzation scheme because Svi/vi = SV2/V2) photon does not mix with 

3) Scalar and pseudoscalar Higgs self energies. We write them as: 

-i^Uip% -i^%{p') (A.25) 

We then have: 

\2 



7=1 



2 -A z]iZ^i (2 ao(mfJ &oK,mfj) 

+ 6 ^ Z^iz'^ (2 aoK) Hml^mD) 
+ 6 ^ ZfZ^- (2 aoKj &oK„m^j) 

3 6 



1=1 



n=l 



6 



m,n=l 



fcJ=l 



fe,«=l 



+ 



\4- E ^M^^ [(2 + 2 m^o - Ml) bo{Ml mlo) 



-w' k=i 



(A.26) 



2 ao(M|) +ao(m^o) + 77^2 <^^'' ao(m^+) - 2 ao(M^) 



+ (2/ + 2m^+-M^) bo{M^,ml, 

^ i ^ i 
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k,l=l 
2 



7=1 



2-^^i^4^ (2aoK) + {P" - boimlml)) 
+ 6 ^ (2 aoiml) + (p^ - 4 m^J 6o(m^^, m^j) 

2 

+ 6 ^ Z^^- (2 ao(m^J + (p^ - 4 m^J 6o(m^^, m^j) 
1=1 

6 

+ E ['^sSTl «o(rniJ + 3 Vs^s'dd M^dJ + 3 ^ssuu M^uJ 

n=l 

- E [3V^^[}V^^^:;&o(rr4,m2^) 

m,n=l 

■ ' k,l=l 

2 / + 2 m^o - Mf) 6o(^|, ) 



g2 /t»fc /ijk 



4 c 



w k=l 

2 ao(M|) +ao(mlfO^^) 



2 



A:+2 



^2 E ^M^i^ 2 ao(M^) 



fc=i 

- f 2 p2 + 2 ml+ - M^] bo{M^, 

k \ k / k 

,4 



+ 4- k(M|) + 2 cl, ao{M^)] - ^^^C^rC^r bo{Ml M|) 



-w 



16 



+ 2 4 6o(M^,M^)l - E o ^S55 ^is-^ bo{mlo,mlo) 



k,l=l 



fe=i 



(A.27) 
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4) Finally we give formulae for the scalar Higgs boson tadpoles which are 
needed to compute counterterms. We write them as: 



(A.28) 



and we have: 



1=1 



4 ^ Z'i a,{ml) + 12 ^ Z'i ao(m^J 



+ 12!^Z|aoK: 



1=1 



+ 



E [^sTl «o(miJ + 3 l-^- «J + 3 V^suu «o(m^j] (A.29) 



n=l 



^Likk Rikk\ 



- 2 E (cscc - cfcc) mCkaoimU - E {csnn - c'snn ) mNMm}^^^ 



Likk „Rikk 



(M|) + 2c2^ao(M^: 



fc=l k=l 

+ 772- 

+ E ^#5 «o(m|o) + ao(m^+) + ^ y^|>^p ao(m|o^^) 

k=l 



Appendix B. Vertex formfactors 

In this Appendix we collect formulae for bare vertex formfactors used in 1- 
loop cross section calculations described in Section 4. Renormalized formfac- 
tors can be obtained be adding to the expressions displayed below the coun- 
terterms defined in Sections 2 and 3. Momentum and Lorentz index conven- 
tions arc shown on Figure (...). Wc neglect momenta in the c-functions argu- 
ments. The notation co(m^,m|,m|) is equivalent to the co(p, g, m^, m|, m|) 
and the same for the other q, Cy . 

B.l Z^-Z*^-scalar vertex 

1. Formfactor proportional to g^^". As in Appendix A, it is very useful 
to define auxiliary functions describing contributions given by the fermion 
triangle loops. Lets define for matter fermions: 

///(m^) = — 3 — 2(p^ + m^)co(m^, m^, m^) — 16c24(m^, m^, m^) 

— 2p^(4cii(m^, m^, m^) + 3c2i(m^, m^, m^)) 

— 2g^(2ci2(m^, m^, m^) + 3c22(^ti^, m^, m^)) 

— 2pg(co(m^, m^, m^) + 2cii(m^, m^, m^) 

-|- 4ci2(m^, m^, m^) -h 6c23(m^, m^, m^)) (B-l) 
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= — 1 + 2(p^ + m^)co(m^, m^, m^) — 2p^C2i(m^, m^, m^) 

+ 2pg (co(m^, m^, m^) — 2c23(m^, m^, m^)) (B.2) 



Contribution of charginos and neutralinos are much more complicated due 
to complex structure of their couplings. Using notation for the vector boson- 
fermion-fermion and scalar /pseudoscalar-fermion-fermion vertices defined in 
( [A.l| , |A.2| ) we can define auxiliary function g^'^J-"^ as: 



9u 



+ 



^Zkl Zmk \ I Imi 



iZklZmk 



cf'cT"'){c' 



Imi 
L 



Jmi\ 
) 



^Zkl„Zmk 



^Zkl„Zmk 



Imi I lm,i\ 



+ {2p'^ +pq)cu{ml,ml,,m^) + {q'^ + 2pq)cu{ml,ml^,m^) 

+ p^C2i{ml, m^, m^) + q^C22{ml, n?^, mf) + 2pq C23(m^, m^, mf) 



nil Tie 



^Zkl Zmk 



_l_ ^Zkl Zmk\( I 



Jimi 



J.mi 



( Zkl Zmk 1 Zkl Zmk\ f Imi , Imiy 
[Cj^ Cy -\- Cy Cj^ ) [Ci^ + ^ 



I 2 / 2 2 2\ 



I /22 2\i2/22 2\i2/22 2\ 

+ pq Ci2{mf^,m^,mi) +p C2i{mj^,m^,mi) + q C22{m^,m^,mi) 
+ 2pq C23{ml, m^, mf) + 2c2i{ml, ml^, ) 



rrim TZe 



(Zkl Zmk 



I Zkl Zmk 



Imi 
L 



Jmi\ 
"R ) 



lmi\ 



^ + iP^+ P(l)cii{ml, m^, mf) 

,2 ™2 



,2 

2 ™2n 



+ (g + pq)ci2{m,^,m^, nil) +p C2i{m,^,m^, Till) + q C22{mi^,m^, nil) 



+ 2pq C23{ml, nil^,nif) + 2c2^{ni\, m^, m]) 
Then we have: 



(B.3) 



E 



4,2 „2 
^^W^W K=l 



72i 



^^2 

2 



4 



+ 4 ^ g^i{ml^,ml^,ml;^ 



+ 2 E ^H("^k'"^L'"^?^,! 



k,Lm=l 



k,Lm=l 
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+ 8 Y: [^e(\41^V|-V^^-i)c24(mi^,mi^,mij 



k,l.rn=l 



+ 3 7^e(\/|t/^\/^^■™^\/-t) c,,{ml^,ml^,mlj 
+ 3 TZeiV^'^^Vl^^V^^'^) C2^{ml^,ml^,mlJ 



o4 2 



2 2 ^ 



g2 2 
g4 2 

"^SwCw j,k,l=l ' ''+^ 



+ 



2p2('r.2 _ „2 \2 2 



Sw^'w j=i 



9^2 



e 



4 



+ 7r^C^(bo{p\M^,M^) + bo{{p + qy,M^,M^) 



2c 



2 
W 



2s2 



2g'c22(M^, M^, M^) - 4pg C23(M^, M^, M^) - 4c24(M^, M^, M^)) 
E co(ml,o, M|, M|) 



4ci4 ^4 

+ E C'/e^M C24(m^o,m2o M|) 

^2„ /i\#2 ji#2 ji#2 \ I /'„2 I o„„\„ /'S,f2 !\#2 j\#2 



+ 2p'c2i{M^.r. M,i., Mf^„) + (g^ + 2pq)c,2{Mty, M^y, M^) 
+ 2g2c22(A4,il4,,il4)+4pgc23(il4,il4,il4) 
+ 17c24(M2,, M^,, M^.)) + 36o(g', M^,, M^,) - 1) 
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k,l=l 



+ 



oi'w'-W K=l 



+ 



+ 



+ 



+ 



e [a 



2 \2 2 



'^■^W^W j=l ^ ^ 



rrir 



.9 ^2 E ^sss ^o(?^ m^o , m^o) 

^^W'-W j=l ^ ^ 



B'r E bo{q^ rnlo^^ , m^o^^ ) 



(B.4) 



^'^^W'-W j=l 



" j+2 



2. Formfactor proportional to g'*^''. Lets define: 



mi Tie 



(„Zkl„Zmk 



I Zikl Zmk\f tmi lmi\ 



_l_ / Zfei Zmfe _|_ Zfei Zmfew jmi _|_ Jmi\l r / 2 2 2\ 

+ (,c^ Cy +Cy ji^c^ ; Ci2[m^,m^,mi) 

+ C22{ml,ml^,mf) 



(B.5) 



Then: 



„2 2 ^ 

^W'-W K=l 



E 



+ 2c22(mf^,mf^,mf^)) 



3m: 



72i 



t'2 

+ 3m^K::5(^ + ~ \^w?) i(^i2{ml^,ml^,mlj + 2c22{ml^,ml^,mlJ) 



Vl 



+ 8 E 92A^l,,ml^,mly 

k,l,m=l 



4 E 92i{rn\^.m 

k,l,m=l 



+ 4 ^ [7^e(yi^^^y|I^yJ^l) (ci2K^,mi^,my + 2c22(m^ 



fe,i,m=l 



+ 37^e(yi^^^y|^^y^^^,) Mml^,ml^,mlJ + 2c22{ml^,ml^, m^^)) 
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^^W'-W K=l 

+ 5ci2 im%o , M| , m%o^^ ) + 2c22 (m|^o , M| , m^^o^^ ) ) 

o2 2 



+ 0,2 ^2 E "^M^M^sIs (ci2(?7T.^o , m^o,m^o) 
,4 2 



+ 2c22 (m^o , mlo , m^o ) ) 
eHcl - ^2 2 



+ ^FTl l'^i2 , mjj+ , mjj+ ) 

^W^W j = l 3 3 3 

+ 2c22(m^+,m^+,m^+)) 



Ci, (co(M^, M^, M^) + 2c22(M^, M^, M^)) 



(-2 



+ 3c22(M^,M^,M^)) (B.6) 



4 2 



2- 



- ci2(m^o,m^o M|)) 

J K-\-Z 

9e4 2 

+ ^w^^ii (ci2(M^, M^, M^) + 2c22(M^, M^, M^)) 
3. Formfactor proportional to (fp^. Lets define: 

<■/// 2\ 1 / 2 2 2\ , / 2 2 2\ 

}2a\^'^ ) — T^^v^ ,m ,m ) + cii{m ,m ,m ) 

+ Ci2(m^, m^, m^) + 2c23(m^, m^, m^) (B.7) 

fibi'^'^) = -^Co(m^,m^,m^) + ci2(m^,m^m^) 

+ 2c23(m^m^m^) (B.8) 

/ _ m "7?^) [/^^A;/ Zmfc Zkl Zmk\( Imi lmi\ 
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+ Cii(mfc,m„,m 



+ mi TZe 



f„Zkl„Zm.k 



) 



- {0^4""' + 4'^cr'){cr + c'r)] [ci2{ml ml, mf) 

+ 2c23{ml,ml,mf) 

+ m„ 7^e [(cf 'cf™^ + ci'='ci™^)(cf ^ - c^^ 

I ( Zkl Zmk , Zkl Zmk\( Imi , lmi\'\ \ ( 2 2 2\ 

-+- {c^ Cy -tCy Cj^ )\Cj^ -tCj^ ;j ^Cu{mi^,m^,mi ) 

+ ci2{ml,ml,m^) + 2c23{ml,ml,m]) 



(B.9) 



E 



„2 2 



m 



'7\i 

2 ^i? /j-i; 



3m; 



_R 

-"2 



(ATKJ + (1 - 



;^^)V2T«j) 



,2 /^drf 



3<T^ (/2a K J + (1 - -ywYmK.)) 



k,l,m=l 
6 



A;,Z,m=l 



2 
Li 



k,l'm=l 



3 7^e(\4t/c/^l"c/^5^'[/) (ci2(m2,^,m?,^,m?,J + 2c23(m^„m^^^,m^J) 
3 ne{V^'hj,V^^j,V^^'^) {ci2{ml^,ml^,mlj + 2c23{ml^,ml^,m%)) 

d 3 

Br E (ci2(m?^, ^^4^, "4^) + 2c23{ml^,ml^,mU) 



oi'W'-W k,l=l 

Ci2(m^o, M|, m^o ) + 4cii(m^o, M|, m%o ) + 2c23(m|^o, M|, m|^o )) 

fc i+2 fc l-\-2 k l-\-2 

„2 2 

\kj Alj jrikl 



nJl 2 



2c23{mH0 mHf^rnjjo)) 

3 + 2 I k 

„4 2 



^i? E ^M^M^^ (ci2(m^o,ml,o^^,m^o 



2c23(m^o,mio ,m^o )) 
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+ 



e c 



w 



„2 \2 2 



„2 2 



• -I i i ' . - - 



0„2 



6cn(M2„M2„M2,)) 

CIj (2co(M^, M^, M^) + 4cn(M^, M^, M^^ 



p4|'„2 _ 2 



+ 3ci2(M^,M^,M^)+4c23(M^,M^,M^)) 
+ 



. 4 4 I] C'iJ^l/^M (ci2(m^o,m2o M|) + c23(m^o,mL 



k+2 



(B.IO) 



,4^2 



+ 



e c 



2s 



fClj, (2cn(M^,M^,M^) 



3co(M^,M^,M^) 



+ 9ci2(M^, Ml, Ml) + 18c23(M^, M^, M^)) 



B.2 Z*^-scalar-pseudoscalar vertex 

1. Formfactor proportional to the momentum of the pseudoscalar. In or- 
der to distinguish between scalar and pseudoscalar couplings we have added 
additional subscripts S* or P in the symbols denoting scalar/pseudoscalar- 
2 chargino/neutralino vertices. Lets define ("Tm" in g{^j is artificial - we 
are using uniform convention ( |A.2| ) for the S and P vertices, but actually 
pseudoscalar couplings do not contain i): 

+ P^(cii(^^, rn'^, Tn?) + 2c2i(m^, rn?, rn?)) — g^Ci2(m^, rn?, rn?) 
+ 2pqc2,{m\m\m')+2c2,{m\m\m')) (B.ll) 



+ 



Xm 



^Zkm( Jki jmlj 



Jki rnlj^ 



^Zkmf Jki mlj 
■^RS^LP 



Iki rnlj\ 
LS^RP ) 



X 



/2i 2\/2 2 2\ 2/2 2 2\ 

+ 2p^C2i(m^, m], ml) + 2pq C23(m^, , m^) + 2c24(m^, m^^ ml) 



ho{q^,ml,ml)) 

^Zknifjki mlj 



m^milm 



Iki mlj . Jki mlj 
^RS^RP "T ^LS^LP )) 



C 



Zkm( Jki mlj 
\(^RS^RP 



Cu{ml^,ml,ml) 

^Zkmfjki mlj 



+ mkmjnlm 



I (^RS^LP 



+ d 



Iki mlj^ 
LS^RP , 



Ca 



,Zkm( Jki mlj 



RS^LP 



Jki mlj ^ 



Jki mlj^ 
■-LS^RP , 



/ 2 2 2\ I / 2 2 2\ 



- m m,Tm r^^m^ lki mlj jki mij-. zkm( Iki mij 
irimmi-Lin \ Ly \t.RS^RP ^ ^LS^LP ) ^ ^A {'^RS'^RP 

Co(m^, m^ ml) + Cii(m^, m^, ml)) 



Iki mlj\ 



^Zkm( Jki mlj 



Jki mlj 



X 



X 



X 

(B.12) 
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Then: 



rylj r^li rU/2 \ "^^UK ry2j ry2i rUU/2 \ 



2 4 

k,l,m=l k,l,m=l 
k,l,m=l 

+ 3 Xm(yi^^^y/^^yi-^) (co{ml^,m%,mlJ + 2c^,(ml^,m%,mlJ) 

"^^W k,l=l 

k I k 

3 2 

^SwCw k,l,m=l + 



e' 



/c+2 

5 



+ Pf^ii (l + 6c24(M^,m^+,M^) + g2(2co(M^,m^+,M^) 



j j 3 

+ 2c23(M^,m^+,M^)) +pg(6co(M^,m^+,M^) + ci2(M^,m^+,M^) 

+ 7cn (M^ , , M^) + 2c2i (M^ , , M^) + 2c23 (M^ , , M^) ) ) 



+ e'M^sw^yj^^L (cn(M^, m^+, M^,) + 2co(M^, m^+, M^,)) 

i j 

+ ^AZ5^^ (cn(M^,m^+,M^) + co(M^,m^+,M^)) 
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3 2 

T^-T E C'rA%Vs% {cn{Mlmlo,rn]jo) + 2co{Mlmlo,rnlo)) 

'^^W^W k,l=l Ik Ik 



o5 2 



+ ^^Ci E C^^M (cii(ml,o, M|, Ml) - coimlo, M|, M|)) 

^Sw'^W k=l 

+ (J + 26o(g^M^,m^,)+4c24K+,M^,m^.) 

+ 2m^+ cii , , ml-+ ) + (c2i , Mp,/ , ) 
- CQ{m\+,M'^,m\+) -2cii{m\+,M'^,m\+)) 

3 i j J 

+ g^(2ci2(m^+,M^,m^+) + C22(m^+, M^, m^+) 
J J J J 

+ 4c23(m^+, M^, m^+)) + 2pq{2c2iim%+ , M^, m%+) 

j 3 3 3 

+ C23(m^+,M^,m^+) - Co{m\+,M'^,m\+) - Cii(m^+, M^, m^+))) 



E^^i^.^^<(i + 2&o(?^M|,m^o 



+ 4c24(m^o, M|, rn]jo^^ + 2m^oCii(m^o, M|, m\o^^ 
+ (C2i (m^o , ^1 , ) - Co (m|^o , M| , o ) 

- 2cii(m^o,M|,m^o )) + g^(2ci2(m^o, M|, m^o ) 
+ C22 (m^o , M| , m^jY^ ) + 4c23 (^^o , M| , m|o^^ ) ) 

- 2pg (co (m^o , M| , m^.^^ ) + Cn (m^,, , M| , m^o^^ ) 

- 2c2i(ml^o, M|, m^^o^^) - C23(m^o, M|, m^o^J)) 

+ -^A% {2h,{p\ Ml, ml.) + M|, m^o)) (B.13) 

+ 7^ A'i, {2bo{p\ M^, ml+) + bi{p^, M^, m^+)) 

ZSwCw ^ ^ 



o3 2 



E KA^A'il (6o((p+g)^m^o^^,m|o) 



+ 2bi{{p + q)^,mlo m%o)) 

fe+2 i 

2. Formfactor proportional to the momentum of the scalar. Lets define: 

//■^(m^) = —ho{q^,m^,m^)~hi{q^,rn^,rn^) 

+ P^(cii(m^, m^, m^) + 2c23(m^, m^, m^) + cu^rri^, m^, m^)) 
+ 2pg(ci2(m^, m^, m^) + C22(Tri^, m^, m^)) (B-14) 
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^W^W k,l=l 

+ 2bi{{p + qf,mlo m%o)) 

fe+2 I 



B.3 Photon-Z -scalar vertex 

1. Formfactor proportional to g'^". Contribution from the fermion triangle 
loops can be expressed in terms of f[l{rn?)-, ///(^^) and (simpler then in 
the Z° case) auxiliary function gl{. 
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Using these functions one can write down expression for the G\ formfactor 
as: 
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2. Formfactor proportional to q^q". Lets define: 



g'^ = ruk Tie 

+ C22{ml,ml^,ml) 
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B.4 Photon-scalar-pseudoscalar vertex 

1. Formfactor proportional to the momentum of the pseudoscalar. Lets 
define: 
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+ 4c23(m^+, M^,m^+)) + 2pg(2c2i(m^+,M^,m^+) 
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2. Formfactor proportional to the momentum of the scalar. Lets define: 



dlOij 



Xm 



Iki klj , Jki klj 
^RS^LP "T ^LS^RP 



^ : 2 (2 2 2\ 



+ p^C2iiml,m^,ml) + (p^ + 2pq + 2ml)ci2iml,m^,ml) 
+ {q^ + 2pq)c22{ml, mf, ml) + 2(/ + pq)c23{ml, mf, ml) 
+ 4:C24{ml,ml,ml) -bi{q^,ml,ml) 



— mim^Im 



Jki klj 
'-RS'-RP 



+ d 



.Iki klj 
LS'-LP 



(B.25) 



Co{ml, mf, ml) + 2ci2(m^, mf, ml)) 



Then: 



k,l=l 



- 2ej: [Xm(y/^KrL') (coK,,mL,miJ+2ci2K^,mL,miJ) 

fe,m=l 



+ 2 Im{Vi,'j^V^^'^) {co{mlj^,mlj^,mljJ+2c,2{mlj^,mlj^,mlJ) 
+ TmiV^'^^V^tj,) {c^{m\,ml^,m\) + 1c^2{m\,ml^,m\)) 
- e'MwV^''^_Ci2iM^,m%+,M^) 



+ 



^kiVs+_ (co(m^+,m^+,m^+) + 2ci2(m^+,m^+,m^+)) 

SW I— I k I k k I k 



2s 



w 



-A% {l + &C2A{M^,m]j+,M^)+p\&CQ{M^,m]j+,M^) 



+ 7cii(M^, m^+, M^) + Ci2(M^, m^+, M^) + 2c2i(M^, m^+, M^) 

J 3 3 

+ 2c23(M^,m^+,M^)) +M(2co(M^,m^+,M^) + Cn(M^,m^+,M, 

3 3 J 

- ci2(M^, m^+, M^) + 2c22(M^, m^+, M^) + 2c23(M^, m^+, M^))) 



+ 



4s2 



(co(M^, M^, M^) + 2ci2(M^, M^, M^) 
Ai^^^-^ (ci2(M^, m^+, M^) - co(M^, m^+, M^)) 



(B.26) 



+ TT^^M (o - ^^+) + 8c24(m^+, M^, m^+) 

ZS^ Z 3 3 3 

+ 2m^+ ci2(m^+,M^,m^+) +p^(c2i(m^+,M^,m^+) 



47 



- Co(m^+,M^,m^+) - 2ci2(m^+, M^, m^+)) 

+ g^(2ci2(m^+, M^, m\+) + 5c22(m^+, M^, m^+)) 

3 J 3 3 

- 2pq{co{mlj+, M^, m^+) - Cii(m^+, M^, m^+) 

3 3 3 3 

+ 2ci2{ml+,M^,ml+)-3c23{ml+,M^,ml+))+2boiq^,M^,m%+)) 

3 3 3 3 3 

References 

[1] J. Ellis, G. Ridolfi, F. Zwirner, Phys. Lett. 257B (1991), 83 
H.E. Habcr, R. Hcmpfling, Phys. Rev. Lett. 66 (1991), 1815 
Y. Okada, M. Yamaguchi, T. Yanagida, Prog. Theor. Phys. Lett. 85 
(1991), 1; Phys. Lett. 262B (1991), 54 

R. Barbieri, M. Prigeni, M. Caravaglios, Phys. Lett. 258B (1991), 167 

[2] S.R Li, M. Sher, Phys. Lett. 140B (1984), 339 

J.F. Gunion, A. Turski, Phys. Rev. D40 (1989), 2325, 2333; Phys. Rev. 
D39 (1989), 2701 

M. Berger, Phys. Rev. D41 (1990), 225 

[3] J. Ellis, G. Ridolfi, F. Zwirner, Phys. Lett. 262B (1991), 477 

A. Brignole, J. Ellis, G. Ridolfi, F. Zwirner Phys. Lett. 271B (1991), 
123 

[4] R. Barbieri, M. Frigeni Phys. Lett. 258B (1991), 395 

M.A. Diaz, H.E. Habcr, preprint SCIPP-91/14, November 1991 
J. Kalinowski, preprint TUM-T31-14/91 (1991) 

P.H. Chankowski, S. Pokorski, J. Rosiek Phys. Lett. 281B (1992), 100 

[5] P.H. Chankowski, S. Pokorski, J. Rosiek, Phys. Lett. 274B (1992), 191 

[6] P.H. Chankowski, S. Pokorski, J. Rosiek, Phys. Lett. 286B (1992), 307 
P.H. Chankowski, S. Pokorski, J. Rosiek, in Proceedings of the XXVI 
International Conference on High Energy Physics Dallas, August 1992 

[7] A. Brignole, Phys. Lett. 281B (1992), 284 

[8] M. Bohm, W. HoUik, H. Spiesberger Fort. Phys. 34 (1986), 687 

[9] P.H. Chankowski, S. Pokorski, J. Rosiek, preprint MPI-Ph-92/116 

[10] See e.g. A.B. Lahanas and D.V. Nanopoulos, Phys. Rep. 145 (1987), 1 
M. Olechowski and S. Pokorski, Phys. Lett, bf 214B (1988), 393 
G.F. Giudice and F. Ridolfi, Z. Phys. C41 (1988), 447 

[11] See e.g. B. Grz^dkowski, J.F. Gunion Phys. Lett. 243B (1990), 301 



48 



J. Rosiek, Phys. Rev. D41 (1990), 3464 

J.F. Gunion, H.E. Haber, G.L. Kane, S. Dawson, The Higgs Hunter's 
Guide (Addison- Wesley, New York 1990) 

W. Siegel, Phys. Lett. 84B (1979), 193 

D.M. Capper, D.R.T. Jones and P. van Nieuwenhuizen Nucl. Phys. 
B167 (1980), 479 

W. HoUik Fort. Phys. 38, (1990) 165 

J. Kalinowski and H.P. Nilles, Phys. Lett. 255B (1991), 134 

See e.g. A. Djouadi, J. Kalinowski and P.M. Zerwas, Exploring the SUSY 
Higgs Sector at e+e~ Linear Colliders: A Synopsis, preprint DESY 92- 
168, November 1992 

Z. Kunszt, F. Zwirner preprint CERN-TH.6150/91 (1991) 

G.'t Hooft, M. Veltman Nucl. Phys. B153 (1979), 365 
G. Pasarino M. Veltman Nucl. Phys. B160 (1979), 151 
A. Axelrod Nucl. Phys. B209 (1982), 349 



49 



Figure Captions 



Figure 1. Diagrams contributing at the 1-loop level to the processes e^e — > 
Z^hP{H^) and e+e" A^hP{H^), included in the calculations of Sec- 
tion 4. 

Figure 2a. Allowed mass regions for and for the different values 
of iTLt and supersymmetric parameters. Supersymmetric parameters 
(M,q,M,i,Mg„„,/x, A) read (in GeV): for rrit = 120 GeV: solid (1000, 
300, 1000, 500, 0), long-dashed (200, 100, 60, 50, 0), medium-dashed 
(1000, 300, 1000, 100, 1000), short-dashed (1000, 300,1000, 500, 1000); 
for mt = 180 GeV: solid (1000, 300, 500, 250, 500), long-dashed (200, 
100, 100, 100, 0), medium-dashed (1000, 300, 1000, 100, 1000). 

Figure 2b. Allowed mass regions for hP and and contour hues of fixed 
tan/3. 

Figure 3. Cross sections for the hghter scalar production as a function of 
Mh mass for fixed values of Ma- 

Figure 4a. Regions in the (M4, Mh) plane in which at least one of the cross 
sections e+e~ — > or e+e~ — > is larger than assumed (Tq, for 
the top quark mass mt — 140 GeV and for the different values of the 
(7o and CMS energy. 

Figure 4b. Regions in the {Ma, Mh) plane in which at least one of the cross 

sections e~^e~ — > Z^h'^ or e~^e~ — > A^h^ is larger than assumed (Tq, for 
the top quark mass rrit — 180 GeV and for the different values of the 
(To and CMS energy. 

Figure 4c. Regions in the {Ma, Mh) plane in which at least one of the cross 
sections e+e^ — s> Z^H^ or e^e^ A^H^ is larger than assumed ag, for 
the top quark masses rrit = 140 and 180 GeV and for different values 
of the (Jo and CMS energy. 

Figure 5a. Regions in the (Af4,tan/3) plane in which at least one of the 
cross sections e+e" — >• Z^hP , e^e~ — >• APhP or e+e^ — Z'^H'^, e+e~ — > 
A'^H'^ is larger than assumed (To, for the top quark mass rrit — 140 GeV 
and for different values of the (To and CMS energy. 

Figure 5b. Regions in the {Ma, tan f3) plane in which at least one of the 
cross sections e+e" Z°h°, e+e' A°h° or e+e" Z°H°, e+e" 
A^H^ is larger than assumed (To, for the top quark mass mt — 180 GeV 
and for different values of the (To and CMS energy. 



Figure 5c. Regions in the {Ma, tan (3) plane in which at least one of the 
cross sections e+e" Z^hP, e+e" A^hP or e+e" Z^H^, e+e' 
A^H^ is larger than assumed ctq, for the top quark mass mt — 140 GeV 
and rrit — 180 GeV, CMS energy 500 GeV and for different values of 
the (Jq. 
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